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Abstract. It was pointed out bv lPaczvhskil ll987t) that the X-ray luminosity of accreting neutron stars is very sensitive to 
the physical properties of the accretion flow close to the innermost stable circular orbit. The X-ray radiation is dominated by 
that emitted in the boundary layer, where accreted matter hits a star surface. The X-ray luminosity of the boundary layer is 
proportional to the local accretion rate. In this note, we estimate local accretion rate variations from the disk that undergoes 
non-stationary axisymmetric perturbations. The perturbations are given by the poloidal-velocity potential. We obtain a simple 
formula describing the modulation of the accretio n rate for the particular case of global vertical disk oscillations that have 
been recently studied by Abra mowicz et alJl2005l) . 
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1. Introduction 

The quasi-periodic oscillations (QPOs) appear in the light 
curves of more than 20 bright low- mass X ray-b inaries 
(LMXBs) with accreting neutron stars Jvan der Klisl boOO). 
Much attention is attracted to the kilohertz QPOs because 
their frequencies are comparable with orbital frequencies in 
the innermost parts of the accretion disks. The orbital fre- 
quency of a particle orbiting the neutron star of mass M 
at the innermost stable circular orbit (ISC O) is ytscq = 
1580(1 + 0.75j) Hz x IAM Q /M jKluzniaket alJl 990). 

Many models have been proposed to explain the ex- 
citation mechanism of QPOs a nd subsequent m odula- 
tion of the X-ray signal (s ee Ivan der Klisl l2000l and 
iMcClintock & Remilfardl 120031 for a detailed discussion 
of observations and models). Recently, it has been sug- 
gested that the high frequency QPOs arise from a res- 
onance between two oscil lation modes of the innermost 
part of the accretion dis k ( Kluzniak & Abramo wicz 2001; 
Abra mowicz & Kluzniakll200ll) . In the case of the neutron- 
stars sources, the modulation of the X-ray radiation may 
origin ate in the m odulation of the local accretion rate 
( Kluz niak & Abramowiczl2004 . 

In LMXBs that are not pulsars, the magnetic field of the 
neutron star is sufficiently weak, allowing the accretion disk 
to extend down to ISCO. The strongest X-ray radiation then 
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originates in the boundary layer, where accreted material hits 
the star surface. Depending on the star radius R+, the amount 
of energy released in the boundary layer exceeds that radiated 
by the whole disk. It gives about 69% of the total luminosity 
if i?* = 3i? s , or even 86% if i?* = 1.5i? S - 

In this context. Ipaczvnskil i 19871) pointed out that a vari- 
ability of X-ray luminosity of accreting neutron stars may be 
governed by physical properties of the accretion flow close 
to ISCO. In Einstein gravity, the inner edge of the pres- 
sure su pported thick accretion disks is slightly bellow ISCO 
i Abramowiczlfl9 8 5) . The material is accreted from the disk 
through a narrow potential nozzle onto the neutron star. Obvi- 
ously, if the innermost part of the disk is not stationary but is a 
subject to some oscillations then the fine structure of the flow 
at the inner disk edge is significantly changed. This strongly 
affects the accretion rate through the nozzle and the result- 
ing X-ray luminosity of the boundary layer. This scenario is 
in agre ement with the re cent observations of iGilfanov et all 
(2003]) and more recentlv lkevnivtsev & Gilfanovl J2005h that 
strongly point to the fact that neutron-star QPOs are modu- 
lated in the boundary layer. 

In sections [2] and [5] we briefly summarize equations im- 
portant for the disk structure close to ISCO and reproduce 
the calculations of the accretion rate through the inner edge 
of the stationary disk. Then in section|4]we calculate the ac- 
cretion rate from the disk that is subject to nonstationary ax- 
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isymmetric perturbations. We derive a simple formula for the 
accretion rate modulation of a vertically oscillating disk. 

2. Disk structure close to ISCO 

We consider an axisymmetric thick disk made of a perfect 
fluid surrounding a neutron star of mass M. The dynamics of 
the fluid is governed by Euler equation, poloidal component 
of which takes the form 
dv . „ £ 2 Vp 



v ■ V v 



V$ = 0, 



(1) 



dt p 
where the bold-face letters refer to the poloidal part of the 
vectors, a = (a r ,a z ), $ is a gravitational potential, r de- 
notes radial coordinate (we employ the cylindrical coordi- 
nates {r, (j>, z}, with the origin coinciding with the center of 
the star) and p, p and £ are the pressure, density and the an- 
gular momentum of the orbiting flow respectively (in general 
all depend on r and z). The azimuthal component of the Euler 
equation gives conservation of angular momentum, 



— + v ■ W = 0. 

at 



(2) 



We assume that the angular momentum is constant in the 
whole volume of the disk, £(r, z) = to, and that the fluid 
obeys the polytropic equation of state, P = Kp 1+1 / n , where 
K and n are polytropic constant and polytropic index, re- 
spectively. In addition, we assume that the poloidal velocity 
v = (v r ,v z ) can be derived from the potential \- Hence, the 
equation 10 is satisfied automatically and the equation Q 
can be further integrated to Bernoulli equation, 
<9y v 2 

-rff + TT + h + U = const = Us- (3) 
at 2 

Here we introduced the poloidal-velocity potential by v = 
Vx, the enthalpy of the fluid h = f£ dp/p — nKp 1 / 11 and 
the effective potential U = $(r, z) + £q/2t 2 . 

As a model of a strong gravitational field of the star, we 
use the pseudo-Newtonian potential $(r) = —GM/(R — 
Rs), where R = \J r 2 + z 2 and ffa is Schwar zschild radius. 
It was introduced by Paczvriski & Wiita ( 1980) and allows us 
to model general relativistic effects using Newtonian calcula- 
tions with remarkable simplicity. Particularly, it gives a cor- 
rect position of the marginally stable orbit at r — nsco = 
3i?s and well reproduces the Keplerian angular momentum 
of test particles orbiting the star, £k = VGMr 3 / (r — Rg). 
The angular momentum is not a monotonic function of r, as it 
is in Newtonian gravity (Rg — > 0). Instead, it has a minimum 
at the marginally stable orbit. 

The structure of the stationary disk is shown in Figure^ 
The inner edge is at radius r = r- ln , where the angular mo- 
mentum of the flow equals the Keplerian value, £q = £(r- m ). 
The Lagrange point is at coordinates [r[ n , 0]. The equipoten- 
tial surface that crosses itself at the Lagrange point is called 
Roche lobe. The corresponding value of the effective poten- 
tial is denoted Un, The equilibrium configuration exists only 
if the surfa ce of the torus is ins ide the Roche lobe, e.g. when 
U s < Kb.. teoveJl 96.4 lA hramowicz et al.ll 978h 

Otherwise, the dynamical equilibrium is impossible and 
the overflowed matter will be accreted through the potential 
nozzle onto the star. 




Fig. 1. Accretion from a thick stationary accretion disk. The 
position of the disk inner edge and the shape of the equipo- 
tentials are determined by the distribution of the fluid angular 
momentum. Here we consider the simplest case of the con- 
stant distribution, £(r, z) = £q. The plot shows projections of 
the equipotential surfaces to the poloidal plane (solid lines) 
and the distribution of the fluid (shaded region). The matter 
that overflows the Roche lobe (the equipotential surface that 
crosses itself) is accreted onto the neutron star. 

3. Stationary flow 

The station ary accretion rate for Ro che overflow was first cal- 
culated by Koz lowski et alJ dl978t) . who used Einstein's the- 
ory. Here we close ly follow the Newtonian calculations of 
|Xb ramowic We consider a small overflow, so that all 

quantities can be expanded to the second order in the vicinity 
of the Lagrange point L. Particularly, the vertical profile of 
the enthalpy can be expressed as 

d 2 h 

dz 2 



,z) = h* K 

' 2 



V, 



K 2 = - 



(4) 



where h* = h(r i n ,0) denotes a maximal value of the en- 
thalpy on the cylinder r = r m . The linear order does not 
contribute because the flow is symmetric with respect to the 
equatorial plane. The thickness of the inner edge is H = 
\Jlh* I ' k. Cl ose to r = rin the a ccretion flow becomes tran- 
sonic. After lAbramowiczl 111 9851) . we assume that the radial 
velocity of the flow equals to the local sound speed and that 
the vertical component of the velocity is negligible compared 
to the radial one. This significantly simplifies the solution be- 
cause it allows us to express the poloidal velocity using the 
enthalpy, 



(5) 
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The local mass flux through the nozzle is m — pv r = 
pes = h n+1 / 2 /K n (l + n) n n x l 2 and the integration over the 
cylinder r = r m gives the total mass flux in therms of the 
central enthalpy h* 



M= 



2ti 



H 



rh&z 



-H 



=(2tt) 



3/2_ 



,1/2 



1 



K(n + 1) 



r(n + 3/2) (h*) 



★ \ra+l 



T(n + 2) k 



where r(x) is the Euler gamma function. 

In the Bernoulli equation we keep the term v 2 /2 and 
neglect only the time derivative because of stationarity of the 
flow. We obtain 

h+U = (l + 7T- ) h + U = U s . (7) 



2n 

The parameter n that determines the shape of the enthalpy 
profile can be expressed using a derivative of the effective 
potential. This introduces the vertical epicyclic frequency io z 
to the problem. From equation (0 we obtain 

By substituting the equations and l|8} and introducing 
AU = lAn — Us w e finally recover the result obtained by 

lAbramowic"zl lH985). 



M- 



■ A{n) — AU 

L0 Z 



n+l 



1 


n 


i 


_K{n + l)_ 




n+ l/2_ 



1+1/2 



(9) 



(10) 



A(n) ee (2tt) 3 / 2 

r> + 3/2) 

T(n + 2) ■ 

4. A perturbed flow 



Now, we suppose that the disk is disturbed and oscillates. 
In that case, the accretion flow will not be stationary any- 
more and in order to described the flow we must use the 
Bernoulli equation Q in the full form. The presence of the 
'non-stationary' term dx/dt breaks however the correspon- 
dence between the enthalpy and the effective potential. The 
equipotential surfaces and the surfaces of constant enthalpy 
will not coincide anymore. If the oscillations are a small per- 
turbation, we can expand the Bernoulli equation in the vicin- 
ity of the stationary flow considered above. 

We suppose that the velocity potential can be expressed 

as 



X(r,t) = X(o)O) +ex(i)(r,t), 



(11) 



where the subscript '(0)' refers to the stationary flow and 
the dimensionless parameter e characterizes strength of the 
perturbation. We assume e <C 1. Then, using the definition 

v(r, t) = Vx(r, t) we find 



2ev 



2ec. 



d*(i) 
dr 



v (1) + e 2 v 2 {1) 



\ dr 



( dX(i) 
V dz 



(12) 



The enthalpy is also affected by the perturbation. The new 
value can be approximated by an expansion in the parameter 

e 

h = h {0) +eh {1) +e 2 h {2) +O{e 3 ). (13) 

By substituting into the Bernoulli equation and equating 
coefficients of same powers of e, we get 



(6) h {1) = - 



'(2) 



dX(i 



<)t 



1/2 



dr 



dX(X 
dr 

d X{l) 
dz 



(14) 



(15) 



This way all thermodynamic quantities are expressed using 
the poloidal-velocity potential. 

To progress further, we need a particular form of the per- 
turbation X(i){ r it)- This is a difficult global problem that 
often involves numerical calculations. Several authors stud- 
ied it under different simplifications. For example, iBlaej 
i 19851) gives all possible modes (e.g. eigenfrequencies and 
eigenfunctions) of slender-torus oscillations. In this limit the 
size of the stationary torus is small enough that the en- 
thalpy can be approxima ted by a quadratic function in the 
whole torus. Blaes ( 1985) considered th e Newtonian g ravita- 
tional field. Recently, Kluzniak & Abramowicz ( 2002) recon- 
sidered the problem in general relativity potential and pointed 
to the existence of a particular mode when the torus moves 
rigidl y up and down ac ross the equatorial plane (see also 
Abra mowicz et alJl2005l for more detailed calculations). The 
eigenfrequency of this mode is equal to the vertical epicyclic 
frequency lo z . The presence of rigid modes in a torus oscil- 
lations has been found also in recent numerical simulations 
(e.g. lLee et all2004lRubio-Herrera & Leel2005l) . 

In the following, we model vertical disk oscillations by a 
simple ansatz for the poloidal-velocity potential 

X(i) = zv z cosujt, (16) 

where uj is the frequency of the oscillations. Calculating ve- 
locity perturbation, we find 

v (i) = v z e z cos tut . (17) 
Hence, ev z = const can be interpreted as the amplitude of 
the vertical velocity. Equations dl4> and dl5l give 



-v 2 . cos 2 tot. 



hm = zv z lo sin ut, hrfy 
The vertical profile of the enthalpy at r 

h{r- ln , z, t) = h* — k 2 z 2 + ezv z io sin a; 

~v 2 z cos 2 cut + 0(e 3 ) 



reads 



1 ,2„,2 



(18) 



(19) 



that is quadratic in the variable z. The position of the enthalpy 
maximum on the cylinder r — r- m is shifted from z = to 
height Sz(t) given as 

L0V z 



Sz(t) = SZsmujt, 5Z = e — — . 



(20) 



We can interpret SZ as the amplitude of the oscillations. Also 
the value of enthalpy in the maximum differs from the sta- 
tionary case by 



Sh*\ 



; h(r„ 
1 , 



, Sz) — h* 



Sz 2 



:(5Z 2 -6z 2 ) 



0(e 3 ). 



(21) 
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Fig. 2. Left: The vertical profiles of the enthalpy /i(r; n , z) on the cylinder r = r- m during vertical disk oscillations (thin lines). 
The amplitude of oscillations is SZ — 0.2H and we chose the polytropic index of the fluid n = 1.5. The figure captures 
profiles with the enthalpy maxima at Sz = 0, 0.1H and 0.2H. The enthalpy profile for the unperturbed stationary disk is also 
shown (thick line). Right: The modulated accretion rate from the oscillating disk (thin solid line). The accretion rate for the 
stationary disk is plotted by thick line. Time is rescaled by the period of oscillations. For reference we plot also the phase of 
disk oscillations (dotted line) The accretion rate is modulated with twice the frequency of oscillations. 



According to equation the actual accretion rate de- 
pends on the maximal enthalpy as M oc (h*) n+1 . This re- 
lation can be applied also in the case of vertical oscillations 
because the z-dependence of enthalpy on the cylinder r = r- m 
can be approximated by a quadratic function also in this case 
and the oscillations do not contribute to the radial velocity 
of accreted matter. Hence, using equations @, l|8} and (12 U 
and assuming that the frequency of oscillations equals to the 
local vertical epicyclic frequency, u = uj z , we arrive at our 
final result 



SM , . Sh* 

(" + 1 )tt 



Mi 



(0) 



P 



2-2p 



(1+p) 



Hp 



sz 2 
} 1p 



(22) 



where SM = M - M (0) and p = n/(n + 1/2). 

Figure|2]shows the result. The enthalpy profiles h(r- mi z) 
are shown for several values of Sz in the left panel. The am- 
plitude of oscillations is SZ/ H = 0.3. The right panel shows 
the modulation of the accretion rate from the oscillating disk. 
The time is rescaled by the oscillation period, T = 2ir/ui z . 
Finally, the time-averaged accretion rate is given by 



M, 



(0) 



H 2 



(23) 



that is positive for reasonable values of n. 

5. Conclusions 

In this note we studied the accretion rate from a non- 
stationary pressure supported accretion disk that undergoes 
the vertical axisymmetric oscillations. The oscillations were 
modelled by a simple ansatz for the perturbation of poloidal- 
velocity field. We believe, however, that several features 



would be present also in more sophisticated (perhaps numer- 
ical) solutions: (1) the first correction to the stationary accre- 
tion rate is of the quadratic order in both the actual perturba- 
tion Sz and the amplitude SZ. This is probably because of the 
symmetry of the stationary flow with respect to the equatorial 
plane. Hence, the frequency of the modulation must be twice 
the oscillation frequency. (2) The accretion rate is maximal 
when the disk reaches the maximal amplitude Sz = SZ. (3) 
The averaged accretion rate from the periodically perturbed 
flow is greater than the that of the stationary flow. 
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